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Superfluidity is an emergent quantum phe¬ 
nomenon which arises due to strong interactions 
between elementary excitations in liquid helium. 
These excitations have been probed with great 
success using techniques such as neutron and 
light scattering[l, 2]. However measurements 
to-date have been limited, quite generally, to 
average properties of bulk superfluid[l] or the 
driven response far out of thermal equilibrium[3— 
6]. Here, we use cavity optomechanics[7—10] 
to probe the thermodynamics of superfluid ex¬ 
citations in real-time. Furthermore, strong 
light-matter interactions allow both laser cool¬ 
ing and amplification of the thermal motion. 
This provides a new tool to understand and con¬ 
trol the microscopic behaviour of superfluids, 
including phonon-phonon interactions[l], quan¬ 
tised vortices[ll] and two-dimensional quantum 
phenomena such as the Berezinskii-Kosterlitz- 
Thouless transition[12]. The third sound modes 
studied here also offer a pathway towards quan¬ 
tum optomechanics with thin superfluid Aims, in¬ 
cluding femtogram effective masses, high mechan¬ 
ical quality factors, strong phonon-phonon and 
phonon-vortex interactions, and self-assembly 
into complex geometries with sub-nanometre fea¬ 
ture size. 

Elementary excitations, in the form of phonons and 
rotons, are fundamental to both the macroscopic and 
microscopic quantum behaviour of superfluid hehum-4, 
including phenomena such as dissipation[4, 5], quantum 
turbulence [13] and quantum phase transitions[14] Tech¬ 
niques to probe such excitations have been crucial to our 
understanding of superfluids since the 1960s[l|. For in¬ 
stance, neutron and light scattering[l, 2, 15] allow the 
dynamic structure factor to be determined, which quan- 
tihes the dispersion relation, as well as the mean occu¬ 
pancy and correlations. However, such techniques are 
slow compared to the characteristic dissipation rate of 
the excitations, constraining them to average thermo¬ 
dynamical properties of the superfluid and prohibiting 
real-time measurement and control. Real-time measure¬ 
ments have only previously been performed by applying 
an external driving force to individual modes in super¬ 
fluid acoustic resonators [5, 6, 16]. This excites them far 
out of thermal equilibrium and constrains measurements 
to coherent dynamics. 

In cavity optomechanics, the coupling between optical 


and mechanical degrees-of-freedom, and therefore mea¬ 
surement rate, is greatly enhanced by the presence of a 
high quality optical cavity. This has enabled the demon¬ 
stration of a range of quantum behaviour [8, 9, 17]; mea¬ 
surement precision approaching the fundamental limit 
set by quantum uncertainty[10|; and precision sensors 
of mass, acceleration and magnetic flelds[18|. Excita¬ 
tions in superfluids have recently been identified as an 
attractive mechanicai degree-of-freedom[6, 19]; introduc¬ 
ing unique features such as viscosity that approaches 
zero at absoiute zero, quantized rotationai motion and 
vortices[ll|, and strong phonon-phonon interactions[l|.In 
the oniy previous experiment a pressure wave in buik he- 
iium acts as a gram-scaie resonator[6|, with the combina¬ 
tion of high mechanicai quaiity factor and mass provid¬ 
ing a path towards uitra-precise inertiai sensors. How¬ 
ever, the comparativeiy iarge mass presents significant 
chaiienges for the observation or controi of thermai exci¬ 
tations, and the manifestation of quantum effects. 

Here we propose and utiiise an aiternative approach to 
superfiuid optomechanics based on femtogram Aims of su¬ 
perfluid heiium condensed on the surface of a microscaie 
whispering-gaiiery-mode resonator (Fig. la). Super¬ 
fluid Aims form naturaiiy on surfaces due to the com¬ 
bination of uitraiow viscosity and attractive van der 
Waais forces. Excitations in such Aims, known as third 
sound[4, 5, 20, 21], manifest as perturbations to the 
thickness with the restoring force provided by the van 
der Waais interaction. The physicai structure of the 
resonator provides a tempiate for the seif-assembiing 
Aim, acting to confine third sound modes at the mi¬ 
croscaie in two dimensions, with Aim thickness defining 
the third dimension. Optomechanicai coupiing is reai- 
ized via the opticai evanescent fleid[22|, with the Aim 
being naturaiiy iocated in the region of maximum fleid 
strength. Compared with previous third sound experi¬ 
ments that use centimetre-scaie enciosures and capacitive 
measurements, this architecture enabies three orders-of- 
magnitude reduction in mechanicai mode voiume com¬ 
bined with greatiy enhanced readout sensitivity[4, 5, 21|. 

To experimentaiiy reaiize thin-flim superfluid optome¬ 
chanics, a fibre coupied microtoroidai resonator is piaced 
in a iow pressure heiium-4 gas environment within a 
heiium-3 cryostat (Fig. lb). At the pressures used the 
helium gas transitions directly to the superfiuid state at 
approximately 1 kelvin, avoiding the normal fluid phase. 
Van der Waais forces then coat the surface of the sample 
chamber with a film of superfiuid helium. A typical op- 
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FIG. 1. Optomechanics with superfluid helium Aims, (a) Illustration of third sound waves of a superfluid helium film 
coating a microtoroid. The third sound oscillations have phase velocity C 3 = \/3a„duid~^, where a-odm = 2.65 x nm®s“^ 
is the van der Waals coefficient for silica, (b) Shot-noise limited homodyne detection of superfluid helium optomechanics. A 
tapered fibre-coupled microtoroid (major diameter = 70 /im, minor diameter = 7 /rm) is mounted inside the sample chamber 
of a helium-3 cryostat. Helium-4 gas is injected into the sample chamber at a pressure of 69 mTorr (at 2.8 K). AM - amplitude 
modulator, BS - beamsplitter, D - photodetector, NA/SA - network/spectrum analyser. The 90/10 BS probes the optical 
transmission through the microtoroid. At the base temperature of the cryostat the superfluid film Is roughly 10 nm thick. 


tical resonance lineshape is seen at temperatures above 
the gas-to-superfluid transition (Fig. 2a (blue)). Below 
the transition, unstable oscillations appear on the blue 
detuned side of resonance (Fig. 2a (orange, green)), char¬ 
acteristic of optomechanical parametric instability. This 
is observed with as little as 40 nW of optical power, over 
one hundred times lower than would be expected from a 
microtoroid mechanical mode[23]. 

To characterize the mechanical response of the super¬ 
fluid film we lock the laser to an optical resonance of 
the microtoroid. Amplitude modulation allows coherent 
driving of the motion of the superfiuid via photother- 
mal actuation[24]. This arises due to the superfluid 
fountain effect where superfiuids flow towards a localised 
heat source. The modulation is monitored via the phase 
quadrature of the optical field. Sweeping the modula¬ 
tion frequency a dense spectrum of mechanical modes 
is observed (see Supplementary Information), ranging in 
frequencies from 10 kHz to 5 MHz consistent with third 
sound waves confined to length scales on the order of the 
microtoroid dimensions. The modes tend to exist in near 
degenerate pairs, separated in frequency by several hun¬ 
dred hertz. We attribute this to breaking of cylindrical 
symmetry due to scattering by a defect on the micro¬ 
toroid surface. As the cryostat cools the frequency of the 
superfiuid modes is observed to decreeise (Fig. 2b). This 
occurs due to increased condensation of helium into the 
superfiuid film, with thickening of the film resulting in a 
weaker van der Waals mediated restoring force. The low 
temperature plateau in frequency, evident in Fig. 2b, oc¬ 
curs when the majority of gaseous helium in the sample 
chamber has been condensed into the superfiuid phase, 
resulting in a mechanical frequency that can be precisely 
tuned by injecting or evacuating helium gas. 


When coupled to a thermal bath at temperature T, the 
root-mean-square motional amplitude of an oscillator is 
given by equipartition to be Sx = ^/kgT/k where k is 
the spring constant of the oscillator. The motion decorre¬ 
lates over a characteristic timescale of 27r/rm, where F^ 
is the oscillator decay rate. Only if measurement preci¬ 
sion better than Sx is achieved within this time scale is 
it possible to track the thermally driven trajectory of the 
oscillator in phase space. This allows thermodynamical 
fluctuations to be studied and controlled in real-time, and 
the fundamental thermomechanical noise floor of force 
and inertial sensing to be reached. While many tech¬ 
niques have been developed to probe the thermal proper¬ 
ties of superfiuid helium, it has proved difficult to achieve 
this regime. For instance, the light scattering measure¬ 
ments in Ref. [15] averaged photocounts for around 30 
minutes to retrieve the thermal motion spectrum of first 
sound waves, while the recent work of De Lorenzo and 
Schwab [6] remains six orders of magnitude away from 
resolving the thermal motion. 

The reduced mode volume and strong evanescent op¬ 
tomechanical coupling achieved in our architecture com¬ 
bine to greatly enhance the capacity to resolve ther¬ 
modynamical fluctuations. To test whether real-time 
measurements are possible we perform homodyne-based 
phase measurement. Spectral analysis on a high quality 
third sound mode at 482 kHz reveals a thermomechanical 
noise peak characteristic of mechanical oscillations (see 
Figs. 2d-h). It was observed that the mechanical quality 
factor increases substantially with decreasing tempera¬ 
ture (Fig. 2c) consistent with previous observations at¬ 
tributed to phonon-vortex interactions[5, 11], with the 
dissipation rate reaching a minimum of F^/27r = 106 Hz 
at 530 mK. To determine the minimum measurement 

















































3 



Optical detuning (MHz) Temperature (mK) Temperature (mK) 



FIG. 2. Superfluid helium mechanics, (a) Observation of superfluid oscillations as the cryostat cools to base temperature, 
while the laser frequency is scanned and optical resonance tracked. Blue, orange, and green traces are offset vertically and were 
respectively taken at 3 K, 1 K, and 0.6 K. (b) Resonance frequency of a particular superfluid mode versus cryostat temperature, 
with a final frequency of 482 kHz. The solid line is a theoretical fit obtained by modelling the condensation of the helium gas. 
(c) Mechanical quality factor of the superfluid mode versus cryostat temperature. Note above 850 mK the signal to noise ratio 
was too low to accurately measure the quality factor, (d-h) Mechanical spectra of the superfluid mode at 530, 700, 800, 850, 
and 900 mK (from left to right). These measurements were performed with the laser coupled to a microtoroid optical mode at 
A = 1555.1 nm with optical decay rate K/2'7r = 22.3 MHz. 


time required to resolve the thermomechanical motion, 
the signal-to-noise ratio (SNR) of this peak was deter¬ 
mined as a function of measurement duration (Fig. 3a). 
At measurement times greater than ‘I'k(T m the SNR is 
relatively constant, as expected, and reaches a value as 
high as 20.5 dB. When the measurement time is reduced 
below 27^IVm the SNR drops linearly. Extrapolating to 
a SNR of one we And that the superfiuid motion is re¬ 
solvable for measurement times as low as 32 /is; a factor 
of 260 times shorter than lir/Vm (8.2 ms), and there¬ 
fore sufhcient to track the superfluid motion in real time 
(Fig. 3a). To demonstrate this capability, we monitored 
the evolution of the superfluid mode in phase space as a 
function of time from a sequence of measurements each 
having 2 ms duration (Fig. 3b). As can be seen, it was 
possible to track the thermal trajectory of the oscillator 
with precision a factor of 6.2 below the thermomechanical 
noise (Fig. 3c). 

The ability to resolve the thermal motion of the super¬ 
fiuid allows thermometry to be performed on superfluid 
modes. Locking the laser to the cavity resonance, it is 
possible to probe the effect of the optical held on the 
superfluid modes in the absence of dynamical optome¬ 
chanical backaction. Spectral analysis of the modes in 
this regime reveals an increase in mode temperature and 
linewidth with laser power. Thermometry of a micro¬ 


toroid mechanical mode shows that this is not due to bulk 
optical heating (see Supplementary Information). This is 
further evidenced by the scaling of both temperature and 
linewidth with optical power, which was observed to be 
sub-linear. Nonlinear scaling of this form has been ob¬ 
served in other optomechanical systems[25], and can be 
explained by the presence of an independent bath, cou¬ 
pled to the superfluid mode, and driven out of equilib¬ 
rium by the optical held (see Supplementary Information 
for theoretical models and experimental data). 

To experimentally probe the effect of dynamical back- 
action we red-detune the optical held from the cavity res¬ 
onance, observing modifications to the mechanical reso¬ 
nance frequency and dissipation rate (Fig. 4). The com¬ 
bination of radiation pressure and photothermal forces 
provides a mechanism to both cool and heat the su¬ 
perfiuid excitations (Fig. 4a,d). It was found that 
the dynamical backaction is dominated by photother¬ 
mal forces [24]. The strength of these forces is deter¬ 
mined by the spatial overlap between optical and me¬ 
chanical modes. Consequently, mechanical modes of sim¬ 
ilar frequency may experience vastly different photother¬ 
mal forces. This is evident in Fig. 4a where the shaded 
mode at 552.5 kHz has a significantly stronger response 
compared with the adjacent mechanical mode. 

The photothermally induced changes in damping rate 
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FIG. 3. Realtime measurements of superfluid mo¬ 
tion. (a) Signal-to-noise ratio of mechanical spectra of the 
482 kHz mode as the measurement time is decreased. Mea¬ 
surements were made at the base temperature of our cryo¬ 
stat (530 mK). (b) Thermal motion of the superfluid mode in 
phase space. Each orange point is the position of the oscilla¬ 
tor, measured in a time 27r/4rm. The X and Y quadratures 
of the oscillator position are acquired by optimally Wiener 
filtering the homodyne signal, then mixing this down to DC. 
X = W {i (t))) cos (urrit), Y = W{i(t)) sm{uirnt), where i{t) is 
the homodyne signal and W is the Wiener filter. The blue 
circle represents the uncertainty of an individual measure¬ 
ment, defined as the standard deviation of the shot-noise in 
the measurement time. The dark blue line shows an exam¬ 
ple trajectory of the oscillator obtained by making successive 
measurements, tracking its motion in real-time over a period 
of 15 ms. (c) Histogram of the position of the superfluid (or¬ 
ange) and the measurement noise (blue) shown in (b). Both 
the superfiuid position and measurement uncertainty are nor¬ 
mally distributed. Statistics were built up by binning the 
result of 4700 individual measurements taken over 5 s of data 
acquisition. 

and resonance frequency of the 552.5 kHz mode are 
shown in Fig. 4b,c. The linewidth broadens from 
ro/27r = 115 Hz on-resonance to rA/27r = 464 Hz at 
a detuning of A = —0.58/t. This behaviour agrees well 
with photothermal theory with a positive photothermal 
coefhcient (see Supplementary Information), as seen by 


the fit in Fig. 4a-c. From the change in damping rate 
we find that detuning cools the third sound mode by a 
maximum factor of ^ = 0.25. 

J-0 i A 

In contrast to the 552.5 kHz mode, dynamical 
backaction experiments on the 482 kHz mode shows 
linewidth narrowing and spring stiffening with red- 
detuning (Fig. 4d), indicating that here the photothermal 
coefficient is negative, opposing the direction of radiation 
pressure. At a detuning of A = —0.58re the linewidth has 
narrowed to rA/27r = 49 Hz from ro/27r = 137 Hz on- 
resonance. This corresponds to an increase in oscillator 
temperature of ^ = 2.8. Again good agreement with 
theory is achieved as seen by the fit in Fig. 4e,f (see 
Supplementary Information). 

Compared to other cavity optomechanics experiments 
that utilize photothermal coupling, the photothermal ef¬ 
fect observed here is exceptionally fast due to the high 
thermal conductivity of superfiuid hehum[l]. The char¬ 
acteristic time constant can be deduced from the func¬ 
tional form of the photothermal response seen in Fig. 
4b,c and Fig. 4e,f and is found to be approximately 
T fa 600 ns. 

Quasiparticles are believed to play a central role in 
both the microscopic and macroscopic behaviour of su¬ 
perfiuid helium. However, our understanding remains 
incomplete, particularly in two-dimensional films. For 
instance, there remains significant debate about how 
third sound dissipates, with phonon-vortex interactions 
thought to play a crucial role [5]. Phonon-vortex in¬ 
teractions are also important for as-yet untested pre¬ 
dictions such as the formation of Onsager vortices 
as two-dimensional superfluids evolve[26] and the de¬ 
cay of two-dimensional superfiuid turbulence at zero 
temperature[27, 28]. By providing the capacity in thin 
superfiuid films to both resolve and control thermody¬ 
namical motion, and with orders-of-magnitude smaller 
superfiuid volumes than has previously been possible, 
this paper provides a new tool to study problems of this 
kind. 

In general, interactions are enhanced as systems be¬ 
come increasingly confined. In superfiuid helium, for in¬ 
stance, the single phonon-vortex interaction energy scales 
as the inverse-square of the confining length-scale. In our 
devices, this energy is approximately 10 Hz, four orders- 
of-magnitude larger than any previous experiment [4], 
and only a factor of ten smaller than the phonon dis¬ 
sipation rate. As a consequence, our devices offer the 
prospect to resolve vortex-phonon dynamics in real-time 
with as few as ten unpaired vortices. Combined with 
atomically smooth surface finish that reduces the occur¬ 
rence of pinned-vortices, this should allow a new ap¬ 
proach to the study of third sound dissipation mecha¬ 
nisms. 

Third sound dissipation also provides a method to 
probe quantum phase transitions such Eis the Berezinskii- 
Kosterlitz-Thouless (BKT) transition[12]. The BKT 









5 



FIG. 4. Optomechanical heating and cooling of two third sound modes. Detuning the optical field results in cooling 
(a) or heating (d) of the superfiuid excitations. The laser is red-detuned with respect to the cavity for both the 552.5 kHz 
mode and the 482 kHz mode. Whether the mode experiences heating or cooling depends on the relative magnitude of the 
radiation pressure force and the sign of the photothermal constant (/3), which itself depends on the optical and mechanical 
mode, (a) Spectra of two closely spaced mechanical modes around 552.5 kHz with varying optical detuning. The shaded 
mode experiences photothermal induced broadening (/3 > 0), with a mechanical linewidth of Tm.l'i'K = 77 Hz (A/k = 0), 
rm/27r = 334 Hz (A/k = —0.25) and rm/27r = 457 Hz (A/k = —0.5). The difference in photothermal response between 
adjacent modes can be explained by local heating from optical absorption by the defect responsible for the observed mode 
splitting, since the defect can be expected to reside at an antinode of the higher frequency mode, and a node of the lower 
frequency mode, (b-c) Change in mechanical damping rate and mechanical resonance frequency as a function of detuning 
for the 552 kHz mode. A maximum relative frequency shift of 5a;m/27r = —60 Hz is measured at A = —0.60k. Solid line: 
theoretical fit to photothermal broadening, (d) Spectra of two closely spaced mechanical modes around 482 kHz with varying 
optical detuning. The shaded mode experiences photothermal induced narrowing (/3 < 0), with a mechanical linewidth of 
rm/27r = 115 Hz (A/k = 0), Vml^TT = 63 Hz (A/k = —0.25) and rm/27r = 36 Hz (A/k = —0.5) . (e-f) Change in mechanical 
damping rate and mechanical resonance frequency as a function of detuning for the 482 kHz mode. A maximum relative 
frequency shift of 5tJm/27r = 23 Hz is measured at A = —0.60k. Solid line: theoretical fit to photothermal narrowing. Traces 
in a and d are offset for clarity. All detuning measurements were taken with 200 nW of launched optical power. 


transition, while well studied in weakly interacting con¬ 
densates, is not fully understood in dense, strongly in¬ 
teracting systems such a superfiuid helium-4. Phonon- 
vortex interactions play an important role near the phase 
transition[12]. The enhanced confinement and measure¬ 
ment precision reported here will allow BKT to be probed 


with precision beyond that available with current tech¬ 
nologies Furthermore, with a factor of three reduction in 
diameter, it should be possible to reach the strong cou¬ 
pling regime, where the single phonon-vortex interaction 
energy is larger than both phonon and vortex dissipa¬ 
tion rates, with phonons and vortices hybridizing into a 













new form of quasiparticle. This offers the prospect of 
studying a new, previously unexplored regime of BKT 
physics. Such physics would be further augmented by 
the capacity to laser control and cool individual phonon- 
modes demonstrated here. While laser cooling is rou¬ 
tinely achieved in solid- and gas-phase systems such as 
cold atom physics, ion trapping, atomic clocks and op- 
tomechanics, it has not previously been demonstrated in 
superlluid helium, or indeed in any liquid. 

We hnally observe that superlluid hlms may open new 
regimes of cavity optomechanics. Hybridized phonon- 
vortex modes may allow experiments where the optical 
held is strongly coupled to an inherently quantized degree 
of freedom. Furthermore, third sound modes have been 
observed in helium hlms thinner than a single atomic 
layer[21]. This regime is inherently highly nonlinear, with 
mechanical zero point huctuations large compared to the 
him thickness. In sensing applications, atom interferom¬ 
etry with superhuid helium enables precise force and in¬ 
ertial sensing. The ability to resolve the thermomechan¬ 
ical motion of the huid demonstrated here could enable 
superhuid force and inertial sensors that operate at the 
thermal noise limit, a capability which until now has not 
been possible. 


* These authors contributed equally to this work. 
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Two fluid model and third sound 

Mechanical motion of liquid helium is described by the two-fluid model, where the liquid is 
composed of a normal fluid and a superfluid component that intermingle without any viscous 
interaction [1,2]. This model gives several types of motion such as first and second sound. The 
mechanical motion studied here is that of third sound in superfluid helium-4, which refers to 
surface waves occurring on nanoscale films (~ 10 nm thick) [3-6]. In this case the film is so 
thin that the normal fluid component is clamped to the substrate, leaving only the superfluid 
component free to oscillate perpendicular to the surface. Third sound is completely analogous 
to classical waves on water; however, the presence of superfluidity is necessary because the 
waves would be quickly attenuated if it were a normal fluid. 

For third sound waves the restoring force is provided by the van der Waals interaction of 
the helium with the substrate. The resulting mechanical oscillations have a velocity (Cg) that is 
highly dependent on the thickness of the superfluid layer [3]: 



where Ps/p is the ratio of superfluid to total fluid density, and avow is the van der Waals 
coefficient (ctyow = 2.65 x 10^^ nm®s“^ for helium on a silica substrate [7]). 

Theoretical Treatment 

Here we will theoretically explore the situation of generalized optomechanical coupling with 
both radiation pressure and photothermal forces [8-10]. As highlighted in Ref. [10] the pho- 
tothermal effect is not well described by an energy conserving Hamiltonian formalism owing to 
its inherently dissipative nature. This is in stark contrast to the well known dispersive optome¬ 
chanical Hamiltonian, where energy and momentum are conserved by a reversible interaction 
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between optical and mechanical degrees of freedom. 

Using an approach similar to that taken in Ref. [10] the dispersive optomechanical Hamilto¬ 
nian is represented in the Heisenberg picture including coupling to a “bath” of resonators that 
introduces both fluctuations and dissipation into the system. Similarly the photothermal effect 
is introduced as a fluctuating force with temporal correlations generated by the thermal response 
to the random absorption of photons. The quantum Langevin equation for the oscillator position 
x{t) is: 


^eff [x{t) + rmX{t) + Ujl,x{t)_ 


Frp + FpT + Fth 

(3A 
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where A is the absorption coefficient given by the ratio of absorbed to circulating optical power 
and Tt is the thermal response time. The dimensionless parameter /? quantifies the relative 
strength of the photothermal process over the radiation pressure such that Fpp = (3AFpp. As 
will be discussed later the absolute value and even the sign of (3 is strongly dependent on the 
spatial overlap of the optical and mechanical mode. The intracavity field is dispersively coupled 
to the mechanical motion so the quantum Langevin equation for the annihilation operator a{t) 
is : 


a{t) = -[n-i (A° + gx{t))] a(t) + (4) 

where the optomechanical coupling is g, A** is the optical detuning, and the cavity decay is 
K = Kin + Ko- It is useful to uotc that photothermal backaction is enabled by the temporal delay 
of the thermal response Tt in (3), and not the cavity induced delay as in standard dispersive 
optomechanics. As such, for photothermal cooling the condition of being “resolved” considers 
only the thermal response time, hence relaxing the stringent requirements on minimizing the 
optical decay rate. 

To expand the equations of motion into linear and nonlinear components the position and 
intracavity field annihilation operator are expressed as a coherent amplitude with quantum flue- 
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tuations i.e. x = x + 6x{t), a{t) = a + 6a{t) and ain{t) = + 5ain{t). Here we will only 

consider a linear photothermal and radiation pressure interaction so the nonlinear terms in (3) 
and (4) become: 


iga{t)x{t) = ig {ax + 6x{t)a + 6a{t)x) (5) 

hga}{t)a{t) = hg (lap + 2a5X^{t)) (6) 


where 5X^{t) = | (5a’f(t) + (ia(t)) are amplitude fluctuations of the intracavity field. For 
simplicity we have chosen the phase of the intracavity field such that a = a*. The steady state 
equation for the optical field is then: 


a = 


K — iA 


(7) 


where the cavity detuning is also modified by the static displacement, A = A° + gx. The 
linearised equations of motion for the fluctuations are then: 


meff [6x{t) + rmSx{t) + Uj‘^6x{t) 


6a{t) 


6X+{t) + ^ * 6X+{t) 


( 8 ) 


2hga 
+ y^2Tmm^skbT^{t) 

— {k — iA) Sa{t) + iga6x{t) + ^/2K^5ain{t) (9) 


The photothermal integral in (3) has been replaced by a convolution with a Heaviside func¬ 
tion H{t) to preserve causality. Transformation into the Fourier domain (i.e. = 5x) 

yields: 


5x 

6a 


2x{oj)hga6X^ 1 + 
igaSx + \/2Kl6ai 


PA 


1 + iujTt 


( 10 ) 

( 11 ) 


which uses the Fourier identity reducing the functional form of the 

photothermal force to a low pass filter with a comer frequency defined by the characteristic 
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thermalization time. For brevity the mechanical and optical transfer functions have been intro¬ 
duced, x(w)“^ = rueff (cu^ — -h iojVyn) and D{oj) = k + i{uj — A) respectively. From (10) 
it is clear both optomechanical effects arise from intracavity amplitude fluctuations 6X~'~ which 
can be expressed as: 


6X^ = - (5a^(—cu) -h 6a{uj)) 


aSx 


~ ^9o + 0{6ain} 


( 12 ) 

(13) 


XzpF D{uj)D*{-uj) 

where the optomechanical coupling rate has been normalized by the oscillators’ zero point 
motion (xzpf), that is go = gxzpF, and the fluctuations of the injected optical field have been 
grouped into a single term 0{5ain}- In the case of a coherent optical drive this term quantifies 
the additional thermomechanical noise generated by optical vacuum fluctuations, commonly 
known as quantum backaction. Here, we discard the term 0{6ain} since dynamical instabilities 
and thermal noise greatly exceed quantum backaction. Substituting the resulting intracavity 
amplitude fluctuations into (10) gives: 


6x 


x(^ 


—Ag^ I a I ^ujm S xm-eff A 


D{oj)D*{—uj) 

X (cu) \/‘2YjYiTnQffk})T 


1 + 


PA 


1 + ioJTt 




(14) 

(15) 


where we have made the substitutions x|pp = h/2meifUJm- The modified mechanical suscep¬ 
tibility x'(to’) includes the radiation pressure and photothermal optomechanical coupling and 
can be written in the form -h -h iuj\rm + Isolating real 

and imaginary terms gives the modification to the mechanical resonance frequency 5uJm and 
decay rate ^F^ as: 


doJm — 


2Cc^'>77 


1 + 


PA \ 2Koj‘^TtPA 


1 -|- ) 1 -h 


t j 


(16) 


/ 2 I a2 2^ TtPA I o /^1 I 

{k + X — uj ) — ^ + 2k\1 + — ^ 


1 + cuV/ 


(iF™ = -A{uj) 


(17) 
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where A{uj) = is proportional to the magnitude of the cavity response. 

Despite extensive experimental work on linear dispersive optomechanics over the last decade 
the photothermal effect has not received the same level of interest, potentially due to the per¬ 
ceived limitations of its inherently dissipative origin. To illustrate the nature of photothermal 
coupling we take the limit where it dominates over radiation pressure, namely [3A » 1. In 
this situation two distinct regimes can be realised depending on the thermalization rate Tt. As 
previously mentioned the thermal response can be represented by a low pass filter with comer 
frequency 1/rt. If the corner frequency is larger than Um then the dominant effect is a shift in 
the mechanical resonance frequency. Conversely if the comer frequency is less than oom then 
the resulting delay in the photothermal force leads primarily to a modified mechanical decay 
rate. These two scenarios are analogous to the picture of “unresolved” and “resolved” dynamics 
in dispersive optomechanics respectively. 

Experimental Details 

Our experiment is contained inside a sealed sample chamber which is mounted within an 
Oxford Instruments closed cycle helium-3 cryocooler with a base temperature of 300 mK. In¬ 
side the sample chamber the microtoroid is mounted on Attocube translation stages for linear 
positioning relative to the tapered optical fibre, which itself is mounted on a custom made glass 
taper-holder to match the thermal contraction/expansion during thermal cycling. Imaging of the 
microtoroid and taper is achieved via mounted microscopes directed through a set of windows 
fixed into the bottom of the cryostat. A stainless steel access tube connected to the sample 
chamber is used to introduce helium-4 gas into the sample chamber. 

Vibrations from the continuously running pulse-tube cooler (PTC) are transmitted through 
the cryostat to the sample chamber, preventing stable positioning of the taper-toroid separation. 
To circumvent this issue, on-chip stabilization beams are added using photolithography during 
the fabrication of the microtoroid. With one beam located on either side of the microtoroid the 
taper may “rest” upon the beams, making the relative vibration common-mode and stabilizing 
the separation to a high precision. This makes it possible to indefinitely maintain critical cou- 
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FIG. 1: Superfluid mode spectrum. Network analysis measurement showing the spectrum of third 
sound modes present in our resonator at 900 mK. Each peak in the spectrum corresponds to a different 
third sound mode. 

pling even while the PTC is operational. To minimize light scattering out of the taper mode the 
support beams are only 5 /rm wide and selectively thinned from 2 /rm to less than 500 nm. 

Ultrasensitive readout of the phase fluctuations imprinted by the motion of the superfluid 
film is achieved via homodyne detection using a fibre interferometer. The motion of the third 
sound waves was observed using both spectral and network analysis. Fig. 1 shows a typical 
mode spectrum acquired via network analysis, showing the large range of third sound modes 
present in our system. To lock the relative phase angle between the local oscillator and signal 
we generate a 200 MHz amplitude modulation before the microtoroid. Mixing down the AC 
component of the photocurrent at the modulation frequency then provides a phase dependent 
error signal for the interferometer that is filtered and applied to a piezoelectric fibre stretcher. 
In this configuration the DC component also provides a dispersive error signal for the cavity 
lock which is enhanced by the local oscillator, enabling stable operation even with nano Watts 
of optical power in the signal arm of the interferometer. 
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FIG. 2: On-resonance behavionr. Optical power dependence of the mechanical properties of a super¬ 
fluid mode at 522.5 kHz (orange circles) and another at 482 kHz (blue circles), when the laser is locked 
to the cavity resonance. Measurements are performed with the laser coupled to a microtoroid optical 
mode at 1555.1 nm. (a, b) Measured thermal energy of the superfluid mode versus input laser power. 
The blue line is a power-law fit to the 482 kHz mode with scaling oc P®-®. The orange line is the fit 
to the 522.5 kHz mode with scaling oc P® ®^. (c) Mechanical linewidth (rin/27r) versus optical power. 
The blue line is a power-law fit given by rni/27r = 16.0 x P®-®® -|- 19.3 and the orange line is the fit 
rni/27r = 49.2 X P® _|_ 23 . 4 . (d) Relative shift in mechanical resonance frequency versus laser power. 
This shift to lower frequencies as the laser power is increased is likely due to the density of the superfluid 
component decreasing as the film is heated, decreasing the third sound speed (see Eq. 1)[11]. 

Non-equilibrium batb coupling 

As commented in the main Letter, it is well-known that dynamical backaction due to detuned 
laser driving in a cavity optomechanical system can allow heating or cooling of the mechanical 
motion. In this section, we investigate the effects of the optical probe on the superfluid modes 
with the laser locked to the cavity resonance, in the absence of dynamical backaction. Specifi¬ 
cally, spectral analysis was performed on two superfluid modes with resonance frequencies of 
522.5 kHz and 482 kHz, as the coupled laser power was varied from 7 — 250 nW. This measure¬ 
ment allows modifications of the mechanical resonance frequency, dissipation rate, and mode 
energy to be measured as function of optical power. 
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In these experiments we observe an increase in mode energy, equivalent to heating, of the 
superfluid mode. This increase in mode energy is characterized by a power law scaling of po.eo 
for the 482 kHz mode and po.69 522.5 kHz mode (Figs. 2a, b). This non-integer power 

scaling of mode energy is markedly different to the linear power scaling expected if the tem¬ 
perature increase was due to quantum backaction or bulk heating of the superfluid. Concurrent 
to the increase in mode energy, it was observed that for both modes the mechanical linewidth 
increases with increasing laser power (Fig. 2c) and the mechanical resonance frequency de¬ 
creases with increasing laser power (Fig. 2d). It should also be noted that we have observed 
behaviour that we attribute to boiling off the superfluid layer, that is, for high laser powers the 
superfluid modes rapidly broaden and increase in frequency until they are no longer observable. 
However this occurs at laser powers greater than 5 /xW, over an order of magnitude higher than 
the powers used in experiment. 

To further rule out bulk heating effects, measurements were performed on a mechanical 
mode of the microtoroid itself at 1.35 MHz. In this case the thermal energy of the mechan¬ 
ical mode, dissipation rate, and resonance frequency were all found to be independent of the 
laser power when varied from 7 nW to 250 nW whilst on-resonance. After confirming that 
on-resonance probing does not heat the mechanical modes of the microtoroid, we varied the 
temperature of the cryostat from 10 K to 600 mK with fixed optical power (100 nW). From the 
linear fit to data in Fig. 3 we found that the temperature of the microtoroid closely follows that 
of the cryostat. 

To investigate the effect of coupling to a generalized non-equilibrium bath we first assign an 

to the bath, where Sbb is. the correspond¬ 
ing noise spectral density[12]. Since this non-equilibrium bath is independent to the thermal 
bath (Tth oc kBT), then the associated coupling rate is not given by the intrinsic decay rate (Fq), 
but by the generalized expression Fb = -^ {Sbb i^m) — Sbb [12]. Importantly, we 

note that if the effective bath temperature is negative, namely Sbb{—^) > Sbb{^)^ then the 
coupling rate is also negative. Taking the high quality factor limit gives the final temperature of 


effective temperature Tb (w) = ^ 


In 


Sbb{^) 

Sbb{—^) 
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FIG. 3: Microtoroid thermomentry. The mode energy, proportional to the integrated power spectral 
density (PSD), of a microtoroid mechanical mode at 1.35 MHz as the cryostat temperature is increased. 
The solid line is a linear fit to the experimental data showing there is good thermal anchoring of the 
microtoroid to the cryostat. 


the oscillator due to the presence of the thermal bath and the non-equilibrium bath: 

_TTo + TbTb 

-tpinal — --• (is) 

1 0 “T 1 B 

Depending on how the non-equilibrium bath is excited via the optical field, specifically the 
ratio of Sbb{—^) to S'bb(w), it is possible to obtain linewidth narrowing or broadening in 
concert with increased mode energy. Furthermore, it has been reported in Ref. [13] that the 
bath temperature may depend nontrivially on the coupled laser power, resulting in a final mode 
temperature that permits the non-integer power law scaling behaviour we observe in Fig. 2. 
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